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Heavy Quarkonium Physics beyond the Next-to-Next-to-Leading Order 
of NRQCD * 

Alexander A. Penin ^^ 

'^II. Institut fiir Thcorctischc Physik, Univcrsitat Hamburg, Luruper Chaussee 149, D-22761 Hamburg, 
Germany 

In this talk I briefly review the recent progress in calculation of the high order corrections to the parameters 
of the nonrelativistic heavy quark-antiquark system in the effective theory approach. 



1. INTRODUCTION 

The theoretical study of the nonrelativistic 
heavy quark-antiquark system m and its ap- 
plication to bottomonium Q and toponium [|j 
physics is of special interest because it relies en- 
tirely on the first principles of QCD. The sys- 
tem in principle allows for a perturbative treat- 
ment with the nonperturbative effects being well 
under control and with no crucial model depen- 
dence. This makes the heavy quarkonium to be 
an ideal place to determine the fundamental QCD 
parameters such as the heavy quark mass niq 
and the strong coupling constant as- Recently 
an essential progress has been made in the the- 
oretical investigation of the nonrelativistic heavy 
quark dynamics based on the effective theory ap- 
proach Q. The analytical results for the main 
parameters of the nonrelativistic heavy quark- 
antiquark system are now available up to next- 
to-next-to-leading order (NNLO) in the strong 
coupling constant and the heavy quark velocity 
P [|-|l| . The NNLO corrections have turned out 
to be so sizeable that it appears to be indispens- 
able also to gain control over the next-to-next- 
to-next-to-leading order (N^LO) both in regard 
of phenomenological applications and in order to 
understand the structure and the peculiarities of 
the nonrelativistic expansion. In this talk we re- 
view the first steps in this direction and consider 
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two particular classes of N^LO corrections: (i) the 
retardation effects arising from the emission and 
absorption of dynamical ultrasoft gluons by the 
heavy quarks and (ii) the corrections enhanced 
by a power of ln{l/as) which are not generated 
by the renormalization group (RG) running of a^. 

2. EFFECTIVE THEORY OF NON- 
RELATIVISTIC HEAVY QUARKS 

In this section we briefly outline the effec- 
tive theory approach to the nonrelativistic heavy 
quark dynamics. Let us consider the near thresh- 
old behavior of the heavy quark vacuum polariza- 
tion function 



(qf^q-^ - 9tJ.vq^) ^{q^) 
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where j^ = q^y^q is the heavy quark vector cur- 
rent. Its imaginary part is related to the normal- 
ized cross section of qq production in the photon 
mediated e^e~ annihilation at energy s = q^ 
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where Qq is the fractional charge of quark q. 
The behavior of the vacuum polarization function 
near-below the threshold s = 4to^ determines also 
the masses and leptonic widths of the perturba- 
tive heavy quarkonium bound states. 



Near the threshold the heavy quarks are non- 
relativistic so that one may consider the quark ve- 
locity j3 = ^ \ — Ami/ s as a small parameter. An 

expansion in (3 may be performed directly in the 
Lagrangian of QCD by using the framework of ef- 
fective field theory. In the nonrelativistic problem 
there are four different scales |^^: (i) the hard 
scale (energy and momentum scale like mq); (ii) 
the soft scale (energy and momentum scale like 
PrUq); (iii) the potential scale (energy scales like 
P'^niq while momentum scales like PrUq); and (iv) 
the ultrasoft scale (energy and momentum scale 
like 0^mq). The ultrasoft scale is only relevant for 
gluons. By integrating out the hard scale of QCD 
one arrives at the effective theory of nonrelativis- 
tic QCD (NRQCD) [||. If one also integrates 
out the soft scale and the potential gluons one 
obtains the effective theory of potential NRQCD 
(pNRQCD) which contains potential quarks and 
ultrasoft gluons as active particles p^. The dy- 
namics of the quarks is governed by the effective 
nonrelativistic Schrodinger equation and by their 
interaction with the ultrasoft gluons. To get a 
regular perturbative expansion within pNRQCD 
this interaction should be expanded in multipoles. 
The corrections from harder scales are contained 
in the Wilson coefficients leading to an expan- 
sion in as as well as in the higher-dimensional 
operators of the nonrelativistic Hamiltonian cor- 
responding to an expansion in j3. The nonrela- 
tivistic expansion in as and (3 provides us with 
the following representation of the heavy quark 
vacuum polarization function near threshold 



n{E) 






C{as)G{0,Q,E) 



(4) 



where E = yfs — 2mq is the qq energy counted 
from the threshold, C{as) is the square of the 
hard renormalization coefRcient of the nonrela- 
tivistic vector current and the ellipsis stands for 
the higher-order terms in f3. G'(x, y, E) is the 
nonrelativistic Green function which sums up the 
[as/PY'' terms singular near the threshold. It 
is determined by the Schrodinger equation which 
describes the the propagation of the nonrelativis- 
tic quark-antiquark pair in pNRQCD 
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where Ji is the nonrelativistic Hamiltonian de- 
fined by 



n 



dl 
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V{x)^Vc{x) + ... .(6) 



Here Vc{x) ~ —Cpas/x is the Coulomb poten- 
tial, Cf — 4/3 is the eigenvalue of the quadratic 
Casimir operator of the fundamental representa- 
tion of the color group, x = |x| and the ellipses 
stands for the higher-order terms in as and (3. 
The Green function has the spectral representa- 
tion 
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where tp^ and ijj\^ are the wave functions of the 
qq bound and continuum states respectively. 

Below the threshold the vacuum polarization 
function of a stable heavy quark is determined by 
the bound state parameters. For the leading or- 
der Coulomb (C) Green function the energy levels 
and wave functions at the origin read 



E^ = - 
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where As — asCFmq/2. Note that for the study 
of the bound-state parameters we have (3 k. ag 
so that we are only dealing with one expansion 
parameter. 

The current status of the theoretical research 
can be summarized as follows. The effective non- 
relativistic Hamiltonian is known in the NNLO 
approximation including the two-loop corrections 
to the static potential ||] and the corresponding 
corrections to the Green function have been ob- 
tained in analytical form in 17-13|. The 0{a1) 
contribution to the hard renormalization coeffi- 
cient has been computed in |^. Thus the com- 
plete NNLO analytical expression for the nonrela- 
tivistic heavy quark vacuum polarization function 
is now available. In the next section we present 
the recent results of the calculation of some N'^LO 
order corrections to the heavy quarkonium pa- 
rameters. 



(b) 



Figure 1. (a) Feynman diagram giving rise to the ul- 
trasoft contribution at N'^LO. The single and double 
lines stand for the singlet and octet Green functions, 
respectively, the wavy line represents the ultrasoft- 
gluon propagator in the Coulomb gauge and the ver- 
tices correspond to the chromoelectric dipole interac- 
tion, (b) Feynman diagram giving rise to the leading 
nonperturbative corrections. The crossed wavy lines 
represent the vacuum fluctuations of the gluonic field. 

3. HEAVY QUARKONIUM PARAME- 
TERS BEYOND THE NNLO OF 
NRQCD 

3.1. Retardation effects 

Let us start with a few general remarks con- 
cerning the structure of higher-order corrections 
in pNRQCD. At next-to-leading order (NLO) the 
only source of corrections is the perturbative cor- 
rections to the hard renormaUzation coefficient 
and the static potential, fn NNLO the higher- 
dimensional operators start to contribute. In 
N'^LO the retardation effects which cannot be de- 
scribed by operators of instantaneous interaction 
enter the game. The retardation effects are in- 
duced by the emission and absorption of virtual 
ultrasoft gluons which "feel" the binding effects 
of quark-antiquark interaction p4| , p5[ . This con- 
stitutes a genuinely new feature which is absent 
in NLO and NNLO and is supposed to be the last 
qualitatively new source of the corrections. 

The leading retardation effects which are un- 
der consideration here arise from the chromoelec- 
tric dipole interaction ^^(rq — Tq) • E of heavy 
quarkonium with virtual ultrasoft gluons. The 
corresponding diagram is shown in Fig. la. An 
interesting distinction from the similar QED pro- 
cess is that after emitting the ultrasoft gluon the 
quark-antiquark pair converts into the color octet 
sate with repulsive Coulomb potential Vij{x) = 
[CaI'^ — CF)ois/x where Ca = 3 is the eigenvalue 
of the quadratic Casimir operator of the adjoint 
representation of the color group. As a conse- 



quence the intermediate color octet states belong 
to the continuum part of the spectrum only. 

The result for this ultrasoft contribution di- 
verges in the ultraviolet (UV) region. This di- 
vergence is spurious. It arises in the process of 
scale separation due to the use of pNRQCD per- 
turbation theory at short distances where it is 
inapplicable. We use the dimensional regulariza- 
tion (DR) with rf = 4 — 2e space-time dimensions 
to handle the UV and infrared (IR) divergences 
of the effective theory which are of the form 1/e" 
(n = 1,2,...) as e ^ [||J2^j2|. Compared 
to the "traditional" NRQCD approach endowed 
with an explicit momentum cutoff and a fictitious 
photon mass to regulate the ultraviolet and in- 
frared behavior 0] this scheme has the advantage 
that the contributions from the different scales 
are matched automatically. In the total N^LO 
result the poles in 1/e in the ultrasoft contribu- 
tion are canceled by the infrared poles coming 
from the hard and soft scale corrections. How- 
ever since these corrections are still unknown we 
subtract the divergent part according to the MS 
scheme. This means that the same scheme must 
be used for the calculation of the hard and soft 
scale corrections. As a consequence the partial 
result for the ultrasoft contribution depends on 
the auxiliary ultrasoft scale /x„s which drops out 
in the total result. 

The corresponding corrections to the Coulomb 
energy levels and the wave functions at the origin 
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Lf = 


-81.5379, 


Lr = 


-5.7675, 


if = 


-37.6710, 


Lt- 


0.7340 , 


if- 


-22.4818, 


Lt- 


2.2326. 



Here L^ and L^ are the QCD analogs of the 
famous QED "Bethe logarithms" . They repre- 
sent a pure retardation effect which cannot be 
reduced to the instantaneous interaction contri- 
bution. For n = 1, 2, 3 their numerical values 
are fE3 



(12) 



3.2. Non-RG logarithmic corrections 

The origin of the logarithmic corrections which 
are not generated by the RG running of a^ is 
the presence of several scales in the nonrelativis- 
tic dynamics. The logarithmic integration over a 
loop momentum between different scales yields 
a power of ln(l//3) and in the approximately 
Coulomb system we have f3 oc as- These cor- 
rections can always be associated with the diver- 
gences of the effective theory p^-|30| . By contrast 
the RG logarithms are well known and may be re- 
summed by an appropriate scale choice. 

In the heavy quarkonium the non-RG loga- 
rithms first arise in NNLO corrections to the wave 
functions at the origin. They are generated by the 
following O{0^) operators in the effective nonrel- 
ativistic Hamiltonian: 
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where S represents the spin of the quark- 
antiquark system and {., .} denotes the anticom- 
mutator. The corresponding corrections to the 
wave functions are proportional to the Coulomb 
Green function at the origin which is UV diver- 
gent and in DR is of the following form 
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where the ellipsis stands for the nonlogarithmic 
contribution. The pole term in Eq. ( |l4| ) is can- 
celed by the ©(a^) IR pole of the hard matching 



coefficient C{as). Thus the scale /i in the loga- 
rithm is to be identified with the hard scale niq 
and the sought corrections for the spin one states 
which are of the main interest read [[ll||l3| 
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The residual operators in the NNLO effective 
nonrelativistic Hamiltonian which arc not con- 
tained in Eq. ( |l3[ ) correspond to the purely per- 
turbative corrections to the static Coulomb po- 
tential. The corresponding corrections to the 
wave functions at the origin [^ 13 IJ] contain 
RG logarithms of the form aj. hi"^{p/asmq) (m = 
1,2) which vanish for fi — agniq. This also 
holds in N^^LO for the RG logarithms of the form 
a^ ln™(/x/Q;sTOq) (rn = 1,2,3) because the ultra- 
soft effects enter the stage only in N'^LO so that 
the corresponding running of the strong coupling 
constant at the ultrasoft scale only becomes rele- 
vant in N'*LO. An important point here is that 
starting from NNLO the hard matching coeffi- 
cient C{as) receives a nonvanishing anomalous 
dimension. Therefore starting from N'^LO the 
running of a^ in C{as) should be taken into ac- 
count to match the scale dependence of the wave 
functions at /i = agniq. 

In N'^LO the non-RG leading logarithmic cor- 
rections are produced by the one-loop renormal- 
ization of the O{0^) operators. In dimensional 
regularization the pole part of the correction is 
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where the first three terms contained within the 
parentheses represent the IR divergence while the 
fourth one embodies the UV divergence of the po- 
tential. The IR poles are canceled by the ultrasoft 
contribution with the characteristic scale aj.mq 
and may be read off from |^,Q while the UV 
poles are canceled by the IR poles of the hard 



coefficients and may be extracted from ||3l| , |32| . 
Thus the divergences of the effective theory en- 
dow the operators in the nonrelativistic Hamilto- 
nian with anomalous dimensions. This results in 
the logarithmic corrections to the spin one I = 
energy levels p8|,p9| 
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Note that the IR poles of Eq. (16) introduce a 
factor ln{Ei /Xs) ~ Ina^ while the UV poles con- 
tribute a factor ln(?7ig/As) ~ ln(l/as). 

Some of the operators with singular coefficients 
in Eq. ( p^ ) also result in the contributions to the 
wave functions which is proportional to the sin- 
gular Coulomb Green function at the origin. The 
overlapping logarithmic divergences lead to the 
double pole in 1/e and therefore to the double 
logarithmic contributions which for the spin one 
states read M 
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The calculation of the subleading single log- 
arithmic contributions to the wave functions is 
much more involved and the complete result for 
QCD bound states is not available. However two 
important parts of these corrections are known. 
The Abelian part of the correction to the ground 
sate wave function can be read of the positronium 
lifetime analysis g^ . Not that one has to exclude 
the one photon annihilation contribution from the 
positronium result because this is absent in QCD 
due to the color conservation. The trivial non- 
Abelian contribution originates from the interfer- 
ence between the non- Abelian part of Eq. (15) 



and the 0{as) term in the hard renormalization 
coefficient. The total known single logarithmic 
contribution to the spin one ground state wave 
function reads 
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Figure 2. The normalized cross section R{E) in 
NNLO (solid line) and with the complete leading log- 
arithmic and known part of the subleading logarith- 
mic N'^LO corrections included (dotted line and bold 
solid line respectively). 
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Several attempts have been made to sum up the 
high order RG and non-RG logarithmic correc- 
tions [EgUSS]. The problem however is not com- 
pletely solved yet. 



PHENOMENOLOGICAL 
TIONS 



APPLICA- 



4.1. Top quark phenomenology 

The relatively large electroweak top quark 
width Tt and characteristic scale of the nonrel- 
ativistic Coulomb dynamics aj.mt are consider- 
ably larger than Aqcd and serve as an effec- 
tive infrared cutoff for long distance nonpertur- 
bative strong interaction effects. This makes per- 
turbative QCD applicable for the theoretical de- 
scription of the threshold top quark production. 
At the same time numerically Tt ~ ot^rnt and 
the Coulomb effects are not completely dumped 
by the non-zero top quark width that should be 
properly taken into account. 

In order to analyze the significance of the N^LO 
logarithmic corrections to the cross section R{E) 
we start from the NNLO result of ||l^ and add 
the contributions from Eqs. (|l^,|l|) and Eq. (|9|). 
The results are plotted on Fig. 2. The input 
parameters are taken to be as{Mz) = 0.118, 



mt = 175 GeV and Ft = 1.43 GeV. The soft 
normalization scale of a^ in the nonrelativistic 
Coulomb problem is determined from the condi- 
tion /is = 2as{fJ,s)'mt- The relatively large soft 
normalization point is taken to improve the con- 
vergence on the perturbative series for the spec- 
tral density in NNLO. The instability of the top 
quark is implemented by the complex energy shift 
E -^ E + iVt in Eq. (||). This accounts for 
the leading imaginary electroweak contribution 
to the pNRQCD Hamiltonian. Numerically the 
non-RG logarithmic corrections are comparable 
to the contribution from the N'^LO RG logarithms 
which may be estimated from the renormalization 
scale dependence of the NNLO result. The effect 
of the N'^LO logarithms is twofold. The normal- 
ization of the cross section is reduced by about 
10% around the 15 peak and the energy gap be- 
tween the 15* peak and the nominal threshold is 
decreased by roughly 10%. This partially com- 
pensates the effect of huge NNLO corrections. 

Since the ultrasoft corrections are scheme de- 
pendent we do not include them into the numer- 
ical analysis. However it is quite interesting to 
estimate their numerical importance. If we take 
fius = a'i'mt to cancel the logarithms of as in 
Eqs. (|l^, 11) the ultrasoft correction to the IS* 
peak energy is about —200 MeV and to the IS* 
peak normalization is only about +2% due to 
some cancellations. Note that one power of «<, 
in these equations refers to the ultrasoft gluon in- 
teraction and should be evaluated at the ultrasoft 
scale a^rrit while the two residual powers of ag 
originate from the Coulomb Green function and 
should be evaluated at the soft scale agmt. 

The stability of the perturbation theory for the 
15' peak energy (but not for its normalization!) 
up to NNLO can be manifestly improved by us- 
ing an infrared safe mass parameter instead of the 
pole mass in the analysis p3| , p^ -p0[ . We should 
note that the corrections under consideration are 
nor related to the RG running of a^ and cannot 
be taken into account by the renormalon based 
mass redefinition. On the other hand in some 
cases the behavior of the perturbation perturba- 
tion theory can be improved by using the direct 
relations between the physical observables |]ll|]. 



4.2. Bottom quark phenomenology 

In the case of bottom quark- antiquark produc- 
tion the nonperturbative effects are much more 
significant and one is led to use the sum rule ap- 
proach to get them under control. Specifically 
appealing quark-hadron duality one matches the 
theoretical results for the moments of the spectral 
density 
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For sufhciently large n the moments are saturated 
by the near-threshold region. Then the main 
contribution to the experimental moments comes 
from the T resonances which are measured with 
high precision. On the other hand for n of order 
O (l/ag) the Coulomb effects should be properly 
taken into account on the theoretical side. In or- 
der to analyze the N^LO logarithmic corrections 
to the T sum rules we upgrade the NNLO result 
of |l|] by including Eqs. (0, |l|) and Eq. (|l|). 
We fix the strong coupling constant and focus on 
the determination of the bottom quark mass. At 
present this seems to be the most interesting ap- 
plication of the sum rules. In the bottom-quark 
case Us at the nominal ultrasoft scale seems to 
be too large for a reliable perturbative calcula- 
tion. Thus to be on the safe side we redefine 
the ultrasoft scale to be 2>a1mb so that in both 
cases as = 0.34. We find that the inclusion of 
the N'^LO logarithms in the sum rules leads to a 
reduction of the extracted MS mass value by ap- 
proximately 150 MeV for moderate values of n, 
5 < n < 15. The result essentially depends on 
/is and /i„s because the scale dependence of the 
correction is compensated only by higher-order 
terms. For example, for iig ~ fius ~ "mb the per- 
turbative series for the moments converges much 
better and the correction to the MS mass is only 
about —30 MeV. This result implies that the un- 
sertainty of the value '7ij^(rnj^) w 4.2 GeV 
obtained from NNLO analysis of T sum rules 
|ll|,|l|,^|0) is at least about 100 MeV. 

In the case of the T(15) meson the local duality 
is expected to work though the nonperturbative 
effects are much more important than in the sum 
rules approach. If we assume the T(15) meson to 
be a perturbative system the results of Sect. 3 can 



be applied to compute the corrections to its mass 
which is determined by the binding energy and 
leptonic width which is determined by the wave 
function at the origin. In this way we find that the 
N'^LO logarithmic correction to the T(IS') mass 
is about —170 MeV and to its leptonic width is 
about —70% (the leading and subleading loga- 
rithms are approximately equal). The ultrasoft 
correction to the T(IS') mass is about —110 MeV 
and to its leptonic width is about —10%. 

It is interesting to compare the perturbative ul- 
trasoft contribution to the leading nonperturba- 
tive contribution of the gluonic condensate due 
to vacuum fluctuations of the gluonic field at the 
scale Aqcd which is generated by the similar di- 
agram with the broken gluon propagator shown 
in Fig. lb in order to conclude how "perturba- 
tive" the heavy quarkonium is. In the case of the 
energy level the leading nonperturbative contri- 
bution is given by p3 
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Using the standard literature value (agG^) ~ 
0.06 GeV we have AEi w 60 MeV which is of 
the same scale as the ultrasoft contribution. 

5. CONCLUSION 

In this talk we reviewed the first steps towards 
the N'^LO analysis of the nonrelativistic heavy 
quark-antiquark system. Two special classes of 
N^LO contributions to the key parameters of 
heavy quark-antiquark bound states were consid- 
ered, namely the non-RG corrections enhanced by 
a power of ln(l/as) and the retardation effects. 
They are not related to the RG running of Ug and 
can be considered as typical representatives of the 
N'^LO corrections. 

We gave the numerical estimates of the above 
corrections for the top and bottom quark systems. 
The corrections turn out to be comparable to the 
NNLO ones and reach 10% in magnitude even in 
the case of top where as ~ 1/10. For T(IS') me- 
son they are out of control. This tells us that the 
NRQCD expansion is not a fast convergent series 
for the physical value of the strong coupling con- 
stant. It is highly desirable to complete the cal- 



culation of N'^LO corrections to get more insight 
of the nonrelativistic heavy quark dynamics and 
the structure of the nonrelativistic perturbation 
theory. Although there is no conceptual problem 
on the theoretical side this analysis is extremely 
difficult from the technical point of view and in- 
cludes, for example, the calculation of the three- 
loop hard matching coefficient and the three-loop 
static potential. 

Finally we would like to mention that the sub- 
leading nonperturbative contributions [Q and 
charm quark mass effects pq] which are relevant 
for the bottom quark physics but numerically are 
not so important as the perturbative corrections. 
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